O © Y O GU + EN HL 


e e 
= © 


e نم‎ 
Je» لت‎ 


G4Course *MyNewCourse :: Construct( ) 1 


G4Course «course = new G4Course(); 
course>title(“Geant4 for Beginners. A crash course”); 
course>author(“Hernán Asorey”); 

course—>email( “asoreyhasmail.com” ) ; 
course>description(“a hands-on Geant4 crash course”); 
course>school(“La Conga Physics”); 
course>site(“github.com/asoreyh/geant4-course”); 
course>year(2023); 

course>duration(4xh); 

course—>license( “CCO 1.0 Universal”); 

return course; 
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Disclaimers 


G4Help(); 
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The jargon 


“Much of the problem we have in 
comprehending specialists in any field is 
in understanding their jargon, and not in 
understanding their ideas. This is 
particularly true for medicine” 


۱ ۱ 


(from above) x y 
(up) 4 (down) x 
SAGITTAL TRANSVERSE CORONAL 
(Right-Left) (Top-Bottom) (Front-Back) 


Cranial 


SINISTRAL POSTERIOR ANTERIOR 
(back) (front) 
Caudal DORSAL N VENTRAL 


(toward spine) (toward sternum) 


Transverse or 
Horizontal Plane 


/ Sagittal Plane 
(Midsagittal Plane 
when at midline) 


Frontal or 
Coronal 


5 in 
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The conventions for this course 


This is not a theoretical course. This is a hands-on 
course. So we will work on the natural environments of 
Geant4: an editor (or IDE) and the (Linux) terminal 
(Windows users>some IDEs included their own terminal) 


Geant4 is enterelly written in CH (mandatory), and 
some bash knowledge is always recommended. So, within 
this course, slides are written in english and 
following the highlight conventions for CH or bash 
within these slides. 


There are some python approaches to Geant4 but they 
are out of the scope of this crash course. 
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Conventions for this course 


e Bash conventions: 
# This is a comment in bash 
$ make # this means run at the command in your CLI as user 


# make # this means to run the command in your CLI as root (sudo) 


e Ct conventions: 
/* This is a CH long comment */ 
// This is a CH short comment, and below there is a typical IDE view 
class MyDetectorConstruction : public G4VUserDetectorConstruction 1 


// your class goes here 


J; 
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The conventions for this course 


OS: I personally recommend any updated (222.04), ubuntu 
flavor (ubuntu, mint, xubuntu, kubuntu, ..). 


However, of course you should use whatever OS you feel 
comfortable, even Windows or iOS. You can also use 


virtualization environments for running G4 (later on this 
course). 
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I strongly recommend using VIm or an IDE for programming 


There are so many possible classes, not all of them following 
an standard naming convention, that it could be helpful to 
take the advantages of an IDE. 


You should explore the several available IDEs. At the end, all 
of them will have the functionalities you expected. 


I always used VIm, but recently started using PyCharm, but 
finally I migrated to VS Code. However, VIm works in every HPC 
environment, and PyCharm has the best LaTeX plugin I've found 
(>Texify). 


My personal 
PC list of VScode 
— extensions 
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// A toolkit for the simulation of the 
passage of particles through matter. 


// Its areas of application include high energy, nuclear and 
accelerator physics, as well as studies in medical and space 
science 


G4Download(“geant4.web.cern.ch/"); 
G4Docs(“geant4.web.cern.ch/docs/"); 


G4AppDocs( “geant4-userdoc.web.cern.ch/UsersGuides/ForApplicationDe 
veloper/BackupVersions/V10.7/html/index.html"); 
G4Examples( “Check the ${geant4_examples}/ dir for extra fun"); 


// IMPORTANT NOTICE 
// I will not spend time showing how to install G4 (and root) > 
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G4 Scope { 


/* 


" 
é 
4 


It is integraly programmed in C/C++ and allow to build applications 
including all the aspects of a Monte Carlo simulation process: 


* 
* 


* 


*/ 


the geometry and building materials of the system; 

the fundamental particles involved and all the physics process 
governing particle interactions; 

the tracking of particles in matter and EM fields; 

the medium/detector response to the passage of these particles 
and their by-products; 

the visualization of the detector and particle trajectories; and 

the capture and analysis of simulation data at different levels 
of detail and refinement. 


Liver Study - Dose 
TPS - Pencil Beam 


Allow absorbed dose 
(E, >[E,]=Gy=J kg) 
calculations 


Allow time-evolving geometries 


/* 


1 
2 
3 
/h 
5 // 
6 
/ 
8 
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A comment about versions. 


By the end of 2022 a new major release, G4 11 
was released (current 11.01p01) 


In this course we will use the latest ہ6‎ 10 
version (10.07p04) 


۱ G4 install uas G4 


# while Geant4 have multiple dependencies, 
some of them are required and some others are 
needed for optional features. 


# Check the official installation guide at 
https: eant4-userdoc.web.cern.ch/UsersGuides/InstallationGuide/html 


# I prepared a bash script for installing the required dependencies, 
root and geant4 at ${HOME}/work. (Ubuntu 23.04) 

# Warning 1: it will change your .bashrc file. 

# Warning 2: it will take time (up to several hours). 


$ curl -Lo install-root-geant4.sh 

https: //ww.dropbox.com/s/ej67f1hc88u7/w1a/install-root-geant4.sh?dl=1 
$ chmod 744 install-root-geant4.sh 

$ ./install-root-geant4.sh 
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# Docker 1s a platform designed to help devs 
build, share, and run modern apps. We handle the 
tedious setup, so you can focus on the code. 


# I prepared two Dockerfiles for this course. Follow the instructions 
and download them from (look for them at the utils directory): 


$ git clone https://github.com/asoreyh/geant4-course.git 
# Otherwise, you can pull the docker image from my docker hub: 


$ docker pull asoreyh/root: latest # root version 6.28.04 (2023) 
$ docker pull asoreyh/geant4: latest # G4 version 10.07.04 (2022) 


# Blank installation. Check the docs! 
# There is also a virtual machine built by the Geant4 collaboration 
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# Important note: You will need to follow these 
steps and provide privileged access to this 
docker to be able to run the QT Geant4 
visualization from docker. 


# open a terminal and enable local access to xhost: 
$ xhost +local:root 


and run the docker (if you don't download the docker images it will 
download them): 


$ docker run --privileged -it -e DISPLAY=$DISPLAY -v 
/tmp/.X11-unix:/tmp/.X11-unix asoreyh/geant4:10.07.04 
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About this course Gal 


This is a hands-on course. During this 4 hours we will code from the 
scratch two Geant4-based applications for: 
o calculate the deposited dose in an organ 
o Simulate an EM shielding 
Hopefully: 
o you will reuse these codes for building your own G4 apps 
o this will give you a global view on building G4 apps 
The final version of codes are available at GitHub, but: 
o we will write them from scratch here. Use GitHub codes only as a 
reference in case of troubles 
o try different geometries, materials, particle beams, .. 


o analyze the differences with your own codes (python, ... ) 


NJ 


About this course 
Introduction 


A first contact 
with G4 and its ° 
examples 


The basics 


Sensitivity and efficiency 
volumes 
Visualization and outputs 


Deposited enersy 


e 0 6 6 


HE blocks Apps 


Structure and common 
practices 
e Detector construction 
e Physics lists 


e Generating particles 


e Calculating the absorbed 
dose in a tissue 

e Simulating the effect of 
shielding 


Geant4 philosophy> Toolkit 


No main code, tools for ° 


Toolkit 


building your own app 


Components (interfaces) 


Actions | OF سر‎ 


The building blocks of Q 
E | | Classes 
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Application ° Your implementation 
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Declarations (.hh) 


G4Help(); 


// In complex codes, it is probable you will include the 
same file at different declaration files. 

// In some cases, this could introduce issues due 
multiple definitions. 

// To avoid that, it is a common and highly recommended 
practice to start any declaration as: 


# ifndef FILE HH // (replace FILE by the declaration) 
# define FILE HH 


/* Continue with your declarations here */ 


# endif 
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CMake í a 


# Configuring geant4 apps is performed by using vr 
$ cmake 66 make 88 make install 


# (or ccmake) 

4 Read the docs for the details, but.. 

4 CMake is an open-source, cross-platform family of tools designed to 
build, test and package software. 

4 CMake is used to control the software compilation process using 
simple platform and compiler independent configuration files, and 
generate native makefiles and workspaces that can be used in the 
compiler environment of your choice. 

# CMake improves the usual (and deprecated) 


$ ./configure 88 make 88 make install 


4 method, by providing an easy way to write Makefiles for compiling 
and installing applications 


ل2 
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A geant4 application... BA 
# .. is actually a cmake project. First we need to create a 
folder 

$ mkdir geant4 


and open this folder in VS Code (or your favorite IDE or 
editor). Then, create a new file called: 


$ vim CMakeList.txt 


# We will use it for configuring our project. 


x2» 
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Our CMakeLists.txt 


# Minimum version of cmake required for compiling or project (I recommended at least cmake v 2.8.12, better v 3) 
cmake minimum reguired (VERSION 2.8.12 FATAL ERROR) 


# The name of the project, we will use "Dose" for the G4 dose calculation applications 


project (G4Dose) 


# Let's tell cmake what packages we reguired or need (not the same!) 
# No UI or VIS, uncomment this and comment the other one 

# find package (Geant4 REQUIRED) 

find package (Geant4 REQUIRED ui all vis all) 


# Include G4 libraries. ${string} are environmental variables 
include ($(Geant4 USE FILE) ) 


# Locate sources and headers 
# S(PROJECT SOURCE DIR) is the directory where CMakeLists.txt is located (S(PWD)) 


include directories (S(PROJECT SOURCE DIR) /inc 
$(Geant4 INCLUDE DIR) ) 

file (GLOB headers $(PROJECT SOURCE DIR) /inc/*.hh) 

file (GLOB sources $(PROJECT SOURCE DIR) /src/*.cc) 


# The name of the excecutable we will use, and the files that will be linked to it 
4 dose is the excecutable, dose.cc the main source, etc 


add executable (dose dose.cc ${sources} ${headers} ) 


# The Geant4 libraries that will be linked to our excecutable... 


target link libraries (dose $(Geant4 LIBRARIES] ) 


# Linking the excecutable with our project 


add custom target (G4Dose DEPENDS dose) 
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G4RunManager() 


/* 


*/ 


This object is the “heart” of any G4 application. It is always mandatory and 
should be defined in your main app.cc code (dose.cc in our example) 


It controls the “flow” of the run 


All the interfaces (G4 toolkit) are defined and provided here: 


* G4VUserDetectorConstruction — geometry construction 

* G4VUserPhysicsList — all your physics is here 

* G4VUserActionInitialization « actions 
* G4VUserPrimaryGeneratorAction < primary particles production 
* G4UserRunAction < optionals... 


G4UserSteppingAction, ..‏ ٭ 


* UlManager, VisManager, .. 


G4RunManager() 
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First step: the main code (dose.cc) 


// 0. I/O operations 
#include <iostream> 

// 1. GARunManager class 
#include "G4RunManager.hh" 
// 2. User interface 
#include "G4Ulmanager.hh" 
tinclude "GAUIExecutive.hh" 
۷ 2 ۵۸ 


tinclude "G4VisManager.hh" 


#include "G4VisExecutive.hh" 


ae main (GANE لاه‎ (eee ece) از‎ 
//1. create the G4RunManager object 
G4RunManager *runManager = new G4RunManager (); 
//5. Initialize the runManager 
// runManager->Initialize(); // uncomment to see what happens 
//2. create the user interfas 


G4UlExecutive *ui = new G4UlExecutive (argc, argv); 
G4Ulmanager *UlManager = G4Ulmanager : :GetUlpointer (); 


//3. visualization manager 


G4VisManager *visManager = new G4VisExecutive (); 
visManager ->Initialize (); 

// 4. start the session - and compile to see what happens 
vl->Sesslonsw ente () 2 

return 0; 


CAE. + os 
Scene tree, Help, History 


Scenetree | Help | History 


Sear 


Command 
+ control 

» units 

۲ ۲ 

> particle 

+ geometry 
۲ tracking 

+ event 
, 

» 

, 

» 
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cuts 
run 
random 
process 
gui 

vis 
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GARunManager obje 


G4RunManager *runManager = new Gi 


the runMa 


ae 


runMenager= 


G4UIExecutiv 


Exec 


G4Ulmanager *UlManager = ) 7424 
: lization manag 
G4VisManager *visManager = new Gi 


visManager ->Initialize (); 


the 


ui->SessionStart (); 
return O; 


4 


dose 


wow c» o m 


Useful tips © 


http;//cern.ch/geant4/ 


Tooltips : 


+ Start a new viewer: 
/vis/open/..." 
For example /vis/open OGL" 

* Execute a macro file : 
Ycontrol/execute my macro file" 


Documentatlon : 


* Visualisation publication : 
The Geant4 Visualization System - A Multi-Driver Graphics System 
„Allison, J. et al., International Journal of Modeling, Simulation, and Scientific Computing, Vol. 4, Su, 
http;//www.worldscientific.com/doi/abs/10. 1142/S17939623 13400011 


Getting Help : 


+ If problems arise, try browsing the user forum to see whether or not your problem has already been encountered. 
If it hasn't, you can post It and Geant4 developers will do their best to find a solution. This Is also a good place to 
discuss Geant4 topics In general. https://cern.ch/geant4-forum 

+ Get a look at Geant4 User support pages: http://cern.ch/geant4/support 


(. 1 (2013) 1340001: 


Output 


Session : 
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First step: the 


Hi Oz 
include <iostream> 
/# Me 


include "G4RunManager.hh" 


I/O operations 


G4RunManager class 


User interface 


include "G4UImanager.hh" 


include "G4UIExecutive.hh" 
Visualization 


include "G4VisManager.hh" 


include "G4VisExecutive.hh" 


e ela CE AIRE göçe, clima are) | 


main code (dose.cc) 


Livvmj LULLL LUIYuL uuac 


eyhficaronte:-/[ pbox/projects/geant4/geant4 des /sr 


HK hee e e hehe e e e ehe hehe e ehe hee e e eee 


Geant4 version Name: geant4-10-07-patch-04 (9-September -2022) 
Copyright : Geant4 Collaboration 
References : NIM A 506 (2003), 250-303 
: IEEE-TNS 53 (2006), 270-278 
: NIM A 835 (2016), 186-225 
WWW : http://geant4.org/ 


KORK OK KOK‏ وچ جرج وج هچ ار ار بد وج و ولج و جر OC‏ از بوچ O OC OC OC DC AC AC EC DEC AC AC OC DC‏ وچ KKK‏ ہو و ار او و و زد 


ee ee eee EEEE ------- G4Exception-START -------- EEEE ------- 
//1. create the G4RunManager object *** G4Exception : 3 
G4RunManager *runManager = new G4RunMana issued by : G4RunManager:: InitializeGeometry 
n ۱ G4VUserDetectorConstruction is not defined! 

//5. Initialize the runManager *** Fatal Exception *** core dump *** 

M -Simitiali O; **** Track information is not available at this moment 
runmanager nitialize Ü uncommer aras Step information is not available at this moment 

5 : 

//2. create the user interfase m. |........ EEEE G4Exception-END 0 EEEE سس‎ 
GAUIExecutive *ui = new GAUIExecutive 
G4UImanager *UlManager = G4Ulmanager : : Ge *** Xteption: Aborting execution *** 


//3. visualization manager 
G4VisManager 
visManager ->Initialize (); 
// 4. 
ui > Sessions E AEE 


return 0; 


*visManager = new G4VisExecutive (); 


start the session - and compile to see what happens 


Aborted (core dumped) 
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G4Course *MyNewCourse :: Construct( ) 1 


G4Course «course = new G4Course(); 
course>title(“Geant4 for Beginners. A crash course”); 
course>author(“Hernán Asorey”); 

course—>email( “asoreyhasmail.com” ) ; 
course>description(“a hands-on Geant4 crash course”); 
course>school(“La Conga Physics”); 
course>site(“github.com/asoreyh/geant4-course”); 
course>year(2023); 

course>duration(4xh); 

course—>license( “CCO 1.0 Universal”); 

return course; 
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ude <iostream> za 


G4RunManager class 


ude "G4RunManager.hh" 


User interface pA qu e $ ./dose 
hice 
lude : . 
yi Before to continue we need to define our 
lude 


“volumes”, i.e., where your app detectors and 
volumes will exist and what are they made of? 


D (always 3 volumes, see next) 

G4 

2) ۰ 

- Create your G4VUserDetectorConstruction 
Z (and register it at your runManager) 

GAN 

6417 

//3. visualization manager نس‎ ted (Sore dered) 

G4VisManager *visManager = new G4VisExecutive (); پا‎ 
visManager->Initialize(); 


// 4. start the session - and compile to see what happens 
ul->Sesslonstare (0) y 
return 0; 
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// Materials 


G4NistManager *nist = G4NistManager::Instance(); 4 


// Materials are made of elements, elements are made of isotopes 
G4MaterialsDocs( “https: eant4-userdoc.web.cern.ch/UsersGuides/ForApplicationDe 
veloper/BackupVersions/V10.7/html/GettingStarted/materialDef.html”); 

e G4Isotope «€ name 8 index. Properties atoms (Z,N(ucleons),molar mass) 

e G4Element < name, index 6 symbol. Properties elements (zn ولا‎ À 

e G4Material < name 8 index, macroscopic properties (p, T, p, state) 


ssp 
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// Materials 


// Materials are made of elements, elements are made of isotopes 


G4MaterialsDocs( “https: //geant4-userdoc.web.cern.ch/UsersGuides/ForApplicationDe 
veloper/BackupVersions/V10.7/html/GettingStarted/materialDef.html"); 


e G4Isotope «€ name 8 index. Properties atoms (Z,N(ucleons),molar mass) 
e G4Element < name, index 8 symbol. Properties elements (zn ولا‎ À 
e G4Material < name 8 index, macroscopic properties (p, T, p, state) 


ssp 


// ects create e molecule or 2 


// create the natural isotopes of H 


G4Isotope new G4Isotope("H", 1, 1, 1.*g/mole); 


G4Isotope *D new G4Isotope("D", 1, 2, 2.*g/mole); 


Need to be created inside a 
function 


// create th lement as a mix of isotopes 


G4Element *elH = new G4Element ("Hydrogen", "H", 2); 
elH->AddIsotope (H, 99.985*perCent) ; 
elH->AddIsotope (D, 0.015*perCent); 


11 ereat che molecule as a macarie 


G4Material *matH2 = new G4Material("H2", 0.08375 * kg/m3, 1); 


matH2->AddElement (elH, 2); 
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// Materials 


G4NistManager *nist = G4NistManager::Instance(); o 


// NIST database > >3000 isotopes, 108 elements, >-300 materials 

14 €2 e Ye Gu Gis. — elements 

// G4 AIR, G4 WATER, G4 CALCIUM CARBONATE, ..  < compounds 

// G4 DNA ADENINE, G4 CYTOSINE, . & biochemical compounds 

// G4 KEVLAR, G4 DACRON, .. — industrial materials 

G4NistMaterialsRef( “https: eant4-userdoc.web.cern.ch/UsersGuides/ForApplicatio 
nDeveloper/BackupVersions/V10.7/html/Appendix/materialNames.html2highlight=nist% 


20materials"); 


// HowTo 


// instanciate the NIST manager 


G4NistManager *nist = G4NistManager:: Instance(); 


G4Material *matWater = nist->FindOrBuildMaterial ("G4 WATER"); 


64715767781" ۰۱016 66۰60 626067 = mise—>Himeorzuu 107156 (MEM CON ۵ 
G4Material *matCaCo3 = nist->FindOrBuildMaterial ("G4 CALCIUM CARBONATE") ; 


// als NIST 03 ٤ 


nist->ListMaterials("all");//simple, compound, hep, space, bio, all 


x 
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// 3 types of volumes: solid, logical, physical 


G4GeometryDocs( “https: eant4-userdoc.web.cern.ch/UsersGuides/ForApplicati ev 
eloper/BackupVersions/V10.7/html/Detector/Geometr eometry.html”); 


e Solid ے‎ what is the shape and geometry (half and half!) of your volume? 


CABO AS ONO RAS و‎ sole > Onan, OSM, ٣ 


e Logical «€ what is made of? 


۶۱۱/۵۳۰ tair = nist 701006 (MEM ATR"); 


G4LogicalVolume *logicWorld = new G4LogicalVolume (solidWorld, air, "logicWorld") ; 


e Physical < where the magic (interactions, ... ) occurs 


G4VPhysicalVolume *physWorld = new G4PVPlacement ) 
0, G4ThreeVector(0.,0.,0.), logicWorld, "physWorld", 0, false, 0, true 
); 
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construction.hh/.cc: materials, volumes 


ifndef CONSTRUCTION HH 
define CONSTRUCTION HH 


always!) 

include "G4SystemOfUnits.hh" 

include "G4PhysicalConstants.hh" 

// 2. NIST class, needed for standard materiales 
include "G4NistManager.hh" 

// 3. detector construction class 

include "G4VUserDetectorConstruction.hh" 

// 6. Volumes: physical, logicals and placements 
include "G4VPhysicalVolume.hh" 

include "G4LogicalVolume.hh" 

include "G4PVPlacement.hh" 


// some standard shapes 


include "G4Box.hh" 
// #include "G4Sphere.hh" 
// 4. choose the name of your class 
class MyDetectorConstruction: public G4VUserDetectorConstruction( 
// 5. The class constructor and destructor, they are public 
public: 
MyDetectorConstruction) ; 
~MyDetectorConstruction) 7 
// 6. the constructor of the physical volume 
uy it will construct the physical volume of your system 
G4VPhysicalVolume *Construct(); 
1 
#endif 


// 1. System of Units and Physical Constants (not mandatory, but... 


#include "construction.hh" 
// as always, first the constructor and destructor 
MyDetectorConstruction:MyDetectorConstruction) {} 
MyDetectorConstruction :~MyDetectorConstruction) {} 
// now the construction function 
G4VPhysicalVolume *MyDetectorConstruction :Construct() { 

// 1. your system is build using materials... using NIST: 


G4NistManager *nist = G4NistManager: :Instance() ; 


// usually (but not always, world is made by Air, ie, G4 AIR) 


G4Material *worldMat = nist->FindOrBuildMaterial("G4 AIR"); 
// 2. your world should have a shape. Let's do a box. 


// It always you have three volumes: 


// solid (the shapes), logical (the materials), physicals (the magic) 


// 2.1 solid: name, x/2, y/2, z/2, use the units!!! 
G4Box *solidWorld = new G4Box("solidWorld", 0.5*m, 0.5*m, 0.5*m); 
// 2.2 logical: assing the material: solid, material, name 
G4LogicalVolume *logicWorld = new G4LogicalVolume( 

solidWorld, worldMat, "logicWorld" 
۳ 
/ 2.3 physical: where the magic occurs: rotation, position(x,y,z), 
/ / associated logical volume, name, motherVolume?, bools (negate 
Jul volumes!), numberOfCopies, check for overlaps (always) 
G4VPhysicalVolume *physWorld = new G4PVPlacement 

0, G4ThreeVector(0.,0.,0.), logicWorld, "physWorld", 0, false, 0, 

true 

۳ 
return physWorld; 


Oo ON XX UO + EN با‎ 
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IIS‏ 


And the new dose.cc 


// 0. I/O operations 


#include <iostream> 


// 1. G4RunManager class 


Si 


#include "G4RunManager. 


// 2. User interface It compiles ۱ ° ) 


#include "G4Ulmanager.hh" 


#include "G4UIExecutive.hh" e æ RE تب‎ e inn s s ا‎ s s s تا‎ s s s s s رت تا تا تب تا ا ا‎ ٦ 


۳۸ 3۰ Visualization asoreyh@caronte:~/Dropbox/projects/geant4/geant4-course/codes/src/build$ make 

[ 33%] Building CXX object CMakeFiles/dose.dir/dose.cc.o 
/home/asoreyh/Dropbox/projects/geant4/geant4-course/codes/src/dose.cc: In function 
“int main(G4int, char**)’: 
/home/asoreyh/Dropbox/projects/geant4/geant4-course/codes/src/dose.cc:26:18: 


| 
#include "G4VisManager.hh" l 
| 
| 
l warning: unused variable ‘UIManager’ [-Wunused-variable] 
| 
| 
| 
| 
| 


finclude "G4VisExecutive.hh" 


// once the detector construction is ready, include it 


Tml dl S ele, Pan 
ی۷۷۷‎ ۴ 26 | G4UImanager *UIManager = G4UImanager: :GetUIpointer(); 


int main(G4int argc, char** argv) { | serras 
[ 66%] Building CXX object CMakeFiles/dose.dir/construction.cc.o 
[100%] Linking CXX executable dose 


[100%] Built target dose 


//1. create the G4RunManager object 


G4RunManager *runManager = new G4RunManager() ; 


// once detector is created, then define it 


// but we are still not ready for init . 
LJ 
runManager->SetUserInitializationnew MyDetectorConstructiorm)); B u E W IL I 1 n 0 t WO r k... ۰ / 


//5. Initialize the runManager 


// runManager->Initialize(); // 


//2. create the user interfase 

G4UlExecutive *ui = new G4UlExecutive(argc, argv); 
G4UImanager *UIManager = G4Ulmanager: :GetUlpointer() ۶ 

//3. visualization manager 

G4VisManager *visManager = new G4VisExecutive(); 
visManager->Initialize(); 

// 4. start the session - and compile to see what happens 
ui->SessionStart(); 


return 0; 


And the new dose.cc ¿my 


#include 


#include 


It compiles! :) 


#include "G4UImanager.hh" 


#inc UlExecutive.hh = — — — — — — — — — — — — — 
t4-course/codes/src/build$ make 

#include "G4V4 005 60 

tınclude Gávi 


urse/codes/src/dose.cc: In function 


#include "G4VisExecuti 


urse/codes/src/dose.cc:26:18: 
ariable] 
r::GetUIpointer(); 


#include "construct 


construction.cc.o 


Dk A 


runManager->SetUser 


ui->SessionStart() ; 


return O; 


نم نم یں + CI‏ ہی O ON‏ 
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Physics, this is why we are here... CA 


// Physics process 


Physics processes describe how particles interact with materials. 
e electromagnetic 

hadronic 

transportation 

decay 

optical 

photolepton hadron 

Parameterisation 


G4PhysicsListDocs( “https: //geant4-userdoc.web.cern.ch/UsersGuides/ForApplication 
Developer/BackupVersions/V10.7/html/TrackingAndPhysics/physicsProcess.html"); 


O CON یت ہی‎ + EN HL 


RR PPP 
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PhysicsList, make it simple.. 2ج‎ 


۳۷٠(2 list‏ ۰۸۰0ی 


e Specify and describe all the particles that will be allowed in the app 
e Specify all the physics process assigned to them 
e Provides a flexible way to setup the physics of your app 


// Yes, but 


e You need to know which process are relevant for the energy scales of your 
application 

e Your physics model accuracy will depend on the physics lists you include, 
DUES 
- include only what you need 
Many physics lists overlap each other 


G4PhysicsListDocs("https: //geant4-userdoc.web.cern.ch/UsersGuides/ForApplication 
Developer/BackupVersions/V10.7/html/UserActions/mandatoryActions.html?highlight= 
physicslistitphysics-lists”); 


O CON OM ہج‎ EN در‎ 


. but no so simple. So — Modular PL 2ج‎ 


There are several ready-to-use modular physics lists 

They were constructed by experts, but they are given as it is... 
Many physics lists overlap each other 

Some current lists, all includes HAD, EM, Decays, neutrons 


FTFP BERT — Current G4 default (collider physics) 

FTFP BERT HP < Idem FTFP BERT but n with En > 20 MeV are treated 
separately by the HP neutron models. Reguires G4NDL and RadiactiveDecay 
Shielding — Simulation of deep shielding (includes HP) 

QGSP BERT — former G4 default. Similar but replaced by FTFP BERT 
QGSP BERT HP < Idem 

G40pticalPhysics & Cherenkov and Scintillation (photon process) 
G4EmStandardPhyscis — EM constructors, see the docs 


0 


0 
© 
(©) 
(©) 
O 


G4PhysicsListGuide(“https: eant4-userdoc.web.cern.ch/UsersGuide 
s/PhysicsListGuide/BackupVersions/V10.7/html/index.html"); 


O ON O 0 + EN نم‎ 


در A PR‏ بر با 
ےم ی EIN‏ 


So, Physics > 3 alternative ways 


e Start from scratch > outside the scope of this course (see the docs) 


e Use Physics constructors > let's see for a simple cases 


o Create a class (suggest to create physics.hh and physics.cc) 
class MyPhysicsList : public G4VModularPhysicsList 


o Construct the list 
۷۷۶۵۸۷3291857 sie 3 6702077976867191) 


RegisterPhysics (new G4EmStandardPhysics )(( // EM constructor 
Fe cis se huysa s (av (epar) (P // or eprleal! processes 
} 
o Register the list in the RunManager (dose.cc) 


runManager ->SetUserInitialization (new MyPhysicsList ()); 


e Or, use Physics List factory > highly recommended, directly at the main app 
#include "G4PhysListFactory.hh" // uncoment if you are using PL factory 
comer. ٠٦ ویر ات‎ ۶ 
(7 asel 
G4PhysListFactory physicsListFactory ۶ 


physicsListFactory .SetVerbose (1); 
G4VModularPhysicsList *physicsList = physicsListFactory .GetReferencePhysList (physicsListName ( ر‎ 


runManager ->SetUserInitialization (physicsList (۶ 


CO NOW + 0 NA 


There we go.. physics.hh/.cc 


ifndef CONSTRUCTION HH #include “aysal lala” 

define CONSTRUCTION HH 

// 1. System of Units (not mandatory, but... always!) 

include "G4SystemOfUnits.hh" // 1. Create the constructor 
// 2. NIST class, needed for standard materiales (idem) MyPhysicsList: :MyPhysicsList 0 ( 
include "G4NistManager.hh" 
Piva ass ns construction alhos RegisterPhysics (new G4EmStandardPhysics()); // 


include "G4VUserDetectorConstruction.hh" only use what you need 


// 6. Volumes: physical, logicals and placements 


include 898 +8 770 RegisterPhysics (new G4OpticalPhysics()); 


include "G4LogicalVolume.hh" } 
include "G4PVPlacement.hh" 


// some standard shapes 
include "G4Box.hh" // 2. Create the destructor 


// #include "G4Sphere.hh" MyPhysicsList::-MyPhysicsList() ( 


) 


// 5. The class constructor and destructor, they are public MS. Register in the main app file -> 


// 4. choose the name of your class 


class MyDetectorConstruction : public G4VUserDetectorConstruction { 


00 
MyDetectorConstruction(); 


^MyDetectorConstruction(); 


// 6. the constructor of the physical volume 

7 it will construct the physical volume of your system 
G4VPhysicalVolume *Construct (); 

y; 

#endif 


And the new lines at dose.cc 


cal ب‎ - ME It compiles! :( 


asoreyh@caronte:~/Dropbox/projects/geant4/geant4-course/codes/src/build$ make 
[ 33%] Building CXX object CMakeFiles/dose.dir/dose.cc.o 
/home/asoreyh/Dropbox/projects/geant4/geant4-course/codes/src/dose.cc: In function 
“int main(G4int, char**)’: 
/home/asoreyh/Dropbox/projects/geant4/geant4-course/codes/src/dose.cc:26:18: 
warning: unused variable ‘UIManager’ [-Wunused-variable] 

26 | G4UImanager *UIManager = G4UImanager: :GetUIpointer(); 


| A 


int main(G4int argc, char** argv) ( 


[ 66%] Building CXX object CMakeFiles/dose.dir/construction.cc.o 

[100%] Linking CXX executable dose 

[100%] Built target dose 
asoreyh@caronte:~/Dropbox/projects/geant4/geant4-course/codes/src/build$ ./dose 


runManager->SetUserInitializationnew MyPhysicsList()); 


r 


UIManager->ApplyCommand("/vis/open OGL"); DS ama uma uma ama um å s s s ص‎ = = = um 


UIManager->ApplyCommand("/vis/viewer/set/viewpointVector 1 1 1); An d WO Y ks ۱ ( 1 k n OW , 1 k n OW ) 


İraw e lume dose 


gc Do 
UlManager->ApplyCommand("/vis/drawVolume" ; = 


finclude ۳ 


int main(G4int argc, char** 


runManager->Se 


UlManager->App 
UlManager->App 


UlManager->App 


[ 33%] Building CXX object CMakeFiles/dose.dir/dose.cc.o 


| 
argv) 1 ۱ asoreyh@caronte:~/Dropbox/projects/geant4/geant4-course/codes/src/build$ make 


1 /home/asoreyh/Dropbox/projects/geant4/geant4-course/codes/src/dose.cc: In function 


Ok, we have volumes, materials and physics.. 
We are almost ready > we need actions! 


Create your G4VUserActionInitialization 


and G4VUserPrimaryGeneratorAction 
(and register them at your runManager) 


p.cc:26:18: 


uild$ ./dose 


O O <a O OQ + EN HL 
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Actions, let's the things evolve 2ج‎ 


// We need two interfaces 


e G4VUserActionInitialization is an interface to create and register the 
G4VUserPrimaryGeneratorAction (mandatory) and other user actions 
o Build() «€ function 


e G4VUserPrimaryGeneratorAction is an interface (action!) to describe how the 
primary particles (injection) should be produced 
o GenerateParticles() < function 
o Typically, but not always > G4ParticleGun 


© 


x 
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There we go.. action.hh/.cc 


ifndef ACTION HH 
define ACTION HH 


include "G4VUserActionInitialization.hh" 


include "G4VUserPrimaryGeneratorAction.hh" 
include "G4ParticleGun.hh" 
include "G4ParticleTable.hh" 


include "G4SystemOfUnits.hh" 


class MyActionInitialization: public G4VUserActionInitialization( 
pubiic; 

MyActionInitialization) ; 

~MyActionInitialization) ; 

virtual void Build() const; // our main function 


1۶ 


class MyPrimaryGenerator : public G4VUserPrimaryGeneratorAction{ 
public: 

MyPrimaryGeneraton) 7 

^MyPrimaryGeneratox); 


virtual void GeneratePrimaries(G4Event*) ; 


private: 

G4ParticleGun *fParticleGun; 
u 

#endif 


#include "action.hh" 
MyActionInitialization:MyActionInitialization) {} 
MyActionInitialization:~MyActionInitialization) {} 

void MyActionInitialization:Build() const { 
MyPrimaryGenerator *generator = new MyPrimaryGeneraton) ۶ 
SetUserAction(generator) ۶ 

} 

MyPrimaryGenerator :MyPrimaryGeneraton) { 

fParticleGun = new G4ParticleGun(1) ; 

} 

MyPrimaryGenerator :~MyPrimaryGeneraton) { 

delete fParticleGun; 

} 

void MyPrimaryGenerator :GeneratePrimaries(G4Event *anEvent) { 
G4ParticleTable *particleTable = G4ParticleTable :GetParticleTable(); 


G4String particleName = "proton"; 


G4ParticleDefinition*particle = particleTable->FindParticle(particleNamg ; 


GAThreeVector pos(0., 0., 0.) > 

G4ThreeVector momdir(0., 0., 1.); 

G4double particleKEnergy = 100. * MeV; 
fParticleGun->SetParticlePosition(pos) ; 
fParticleGun >SetParticleMomentumDirectionmomdir) ; 
fParticleGurn->SetParticleEnergyparticleKEnergy ; 
fParticleGum»SetParticleDefinitionparticle); 
fParticleGun >GeneratePrimaryVertexanEvent) ; 


} 


// and finally inform our application to draw the trajectory -> dose.cc 


And the new lines at dose.cc 


nce e ead n e 


#include "action.hh" 


int main(G4int arge, char** argv) | 


nce Ge s abs eated e e 
n 1 hysi i rea OY: To 


runManager->SetUserInitializationnew MyPhysicsList()); 


UIManager-»ApplyCommand("/vis/viewver/set/autorefresh true; 
UIManager->ApplyCommand("/vis/scene/add/trajectories smooth); 


É 1 ect - jethe en eede 


r, 1 


UlManager->ApplyCommand("/vis/scene/endOfEventAction accumulate"; 


And it works! 


asoreyh@inferno:~/Dropbox/projects/geant4/geant4-course/codes/src/build$ make 

[ 20%] Building CXX object CMakeFiles/dose.dir/dose.cc.o 

[ 40%] Building CXX object CMakeFiles/dose.dir/action.cc.o 

[ 60%] Building CXX object CMakeFiles/dose.dir/construction.cc.o 

[ 80%] Building CXX object CMakeFiles/dose.dir/physics.cc.o 

[100%] Linking CXX executable dose 

[100%] Built target dose 
asoreyh@caronte:~/Dropbox/projects/geant4/geant4-course/codes/src/build$ ./dose 


Bom 


Useful tips vlewer-0 (OpenGLs 
Help | History 


Command 
control 


+ random 
» pre 
material 
» gun 
» gui 
vis 


Output 


Session 


N O <a ہی‎ 0 + EN P 
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And the new lines at dose.cc 


// Some useful uiManager commands 


402 


elek “el 


= 


IManager->Appl 
IManager->Appl 
IManager->Appl 


IManager->Appl 
IManager->Appl 
IManager->Appl 
IManager->Appl 


yCommand ("/vis/open 


OGL"); // open the visualization 


positive particle 


= 


IManager->Appl 


a 


IManager->Appl 


|yCommand ("/vis/model 


|yCommand ("/vis/model 


yCommand ("/vis/viewver/set/autorefresh true"); 


yCommand ("/vis/viewer/set/viewpointVector 1 1 1"); 


yCommand("/vis/ogl/set/displayListLimit 50000"); 


yCommand("/vis/drawVolume"); // draw the volumes 


yCommand("/vis/scene/add/trajectories smooth"); 


yCommand("/vis/scene/endOfEventAction accumulate 300"); 


neutral particles 


۸۲۸ set valew ponte 


7۴7 atra co loe ۱1809117-22 


// draw particle trajectories 


// use smooth traj (magnetic field!) 


// accumulate events display 


negative particles 


ling/trajectories/create/draw 


21 


ling/trajectories/create/draw 


ByParticl 


SID"), // mes Id, me charge 
eID-0/set e- blue"); 


O © <a © یت‎ + EN P= 
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Computational physics 1s beautiful 


NE 
// on 


fincl 


int mi 


Ld$ make 


۱08 ./dose 


fun (water) EN 


Increase 


density, 


Increase 


۳ (GN) (oo) cap Un We) Ips (ce) (on 


° 
٦ 


A 
¬ 


cN 
لی‎ 


mM 
اس‎ 


a 
w 


CON DUFF EN HL 


FTFP BERT HP, 


ye Naa 
// once the action is r 


#include "action.hh" 


int main(G4int argc, ch: 


7 ۷۹ 


// once the physics is 


runManager->SetUserInit 


+8 7 

// after actions ... dr 
UIManager->ApplyCommand 
UIManager-»ApplyCommand 
// for viewing many tre 


UIManager-»ApplyCommand 


GE 71 
) 


n, 


1 MeV in water 


photons 
electrons 


neutrons 


jeant4/geant4-course/codes/src/build$ make 
's/dose.dir/dose.cc.o 
's/dose.dir/action.cc.o 
's/dose.dir/construction.cc.o 
's/dose.dir/physics.cc.o 


jeant4/geant4-course/codes/src/build$ ./dose 


O O <a ON + EN HL 


pop 
h e 


pa‏ با 
O‏ + 


Walt! Up to now: 


We have the basis of an app for simulate many physics 


applications 
e Your main app base code (dose.cc) including: 

o The visManager 

o The uiManager 

o The runManager 
m MyDetectorConstruction 
m MyPhysicsList 
m MyPrimaryGeneratorAction 


m MyActionInitialization 


O O — ON + EN HL 


ES کا‎ 
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Wait! Before changing all your codes 


Git is a wonderful tool. If you don't know how to use it please please 


please learn the basis (see e.g. this official tutorial): 


G4GitTutorial(“https://docs.github.com/en/get-started/quickstart/hello-world”); 


I recommend you to create a new branch and work directly on it. So, in 
your repository, checkout the master branch: 

$ git checkout -b testing 

Here, 'testing' is the name I selected for the branch. You can use 
whatever you prefer. Let's check if we are in the correct branch 

$ git branch 

master 


* testing 


For going back to the master branch, just: 
$ git checkout master <=> 


O CON ی ہی‎ + EN HL 
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Now we are safe, let's code: 2ج‎ 


e Before to continue... you MUST play a lot with the different 


options. You should try, at least: 


O 


O 


0 


Different 
Different 
Different 
Different 
Different 


materials (build your own and/or use the NIST DB) 
particles (check the docs) 

energies (see what happens) 

shapes (check the docs) 


physics lists 


Oo ON یت‎ GU + EN با‎ 


S م‎ 
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Now we are safe, let’s code: 


e Perhaps you introduce a lots of changes in your ‘testing’ branch 
e You can merge them into the master branch or continue developing 
in your testing branch 
e Now we let’s start working in the ‘final’ app. Somethings needs 
to be done. 
o Copy the ‘base’ directory to the new ‘final’ directory 
$ cd /path/to/apps/ 


$ cp -r base final 

o Simple but wrong way to avoid building issues: change the 
CMakelist name in the base app. 
$ cd base 


$ mv CMakeLists.txt CMakeLists.off # use the name you want. <j | 


O ON O 0 + EN نم‎ 
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nm‏ بر با 
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Let's define some bio materials and e 


We want to construct a new volume (box?) of skeletal muscle 


cilarerlal “muscle = mise—>rimelo run ۱07۳ پ00‎ ٢ 5 7٥7 


And add the corresponding solid, logic and physical volume: 


o The 'solid' volume, shapes 


G4Box *solidMuscle = new G4Box("solidMuscle", 0.025*m, 0.025*m, 0.025*m); 


o SO IC CR 
— دا‎ 


o The 'logic' volume, materlâTs'—.., 


G4LogicalVolume *logicMuscle = new G4LogicalVolume (solidMuscle, muscle, "logicMuscle"); 


o The ‘physics’ volume, placement and mother volumes 


This is the ‘mother’ volume 


G4VPhysicalVolume *physMuscle = n G4PVPlacement ) © 


0, G4ThreeVector(0., 0., 0.10") >> logicMuscle, "physMuscle", تس‎ s IN ere); 
Let's use a beam of 10 MeV e (action.cc) 


i" —gmw gw 
G4String particleName - "e-"; 114 mM, ۱٦ 


G4double particleKEnergy = 10. * Mev; N "9 ۹ ۰ E ym 


K 0 


=> 
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A comment about doses 


Absorbed dose > Physical magnitude, the energy deposited in matter 


by ionizing radiation per unit mass 


D= ES 


m 


unit: gray (Gy) > [D] = Gy = J 8 (1 rad= 10 erg/g = 107 Gy) 
Equivalent dose — stochastic health effects by ionizing radiation 


R (biological effectiveness, depends on type and E) in a tissue T, 


Total equivalent dose in — > Absorbed dose by the 
=== PESE -Hr == R WRÜRT- sacs — 


the tissue T tissue T due to exposure 


4 
Sum over all the d I 
Radiation weighting 


7 
radiation types and 


energies involved factor (regulations) 


Unit: sievert (Sv) > [H] = Sv = J kg? (1 rem = 107 Sv) 


A 


comment about doses 2 
SA 


Effective dose > stochastic health risk to the whole body due to 
radiation exposure. It takes the nature and biological response of 


each tissue 


Tissue weighting factor 


Sum over all tissues (regulations) 
> ۱ Equivalent dose absorbed 


> by tissue T 
: — es > — mif 
Effective dose E = ۱ T wr] Hr 


Mass-averaged absorbed 


Eo. 
zd dose 


E = b WT > RURDRT POR QUE 
HACERLO FACIL 


Unit: sievert (Sv) > [H] = Sv = J kg! (1 rem = 102 Sv) | «s 


Weighting factors w, and w, 


Wrixon (2008), 
doi:10.1088/0952-4746/28/2/R02 


Table 2. Recommended radiation weighting factors. 


Radiation type Radiation weighting factor, wg 
Photons 1 
Electrons and muons 1 
Protons and charged pions 2 


Alpha particles, fission fragments, heavy ions 20 


Neutrons A continuous function of neutron energy 


2.5 + 18.2exp[— 
Wn = 4 50417. 0 exp|— 
2.5 + 3.25 exp|- 


In(E, | E, < 7۷ 
In(2E, )?] 1 > E, > ۷ 
In(0.04E,,)?] E, > ۷ 


Radiation Weighting Factors for Neutrons 


=ICRP 60 table 
ICRP 60 equation 
*"|CRP 103 
US NRC 


5 
Š 
q 
u 
o 
5 
E 
> 
š 


0 
.000001 .0001 01 1 100 
Neutron Kinetic Energy (MeV) 


w AS Still controversial, ce 98 06 77 


Harrison (2021), 


doi:10.1088/1361-6498/abe548 


Table 1. Summary of ICRP publication 103 nominal cancer risks and detriment for uni- 
form whole-body exposure to gamma rays for the whole population, 0-84 years of age 
(from table A.4.1, publication 103, Annex A). 


Nominal risk 


coefficient 

(cases per 10000 Relative Tissue weighting 
Tissue persons per Gy) Detriment detriment* factor, wr 
Oesophagus 15 13.1 0.023 0.04 
Stomach 79 67.7 0.118 0.12 
Colon 65 47.9 0.083 0.12 
Liver 30 26.6 0.046 0.04 
Lung 114 90.3 0.157 0.12 
Bone surface 7 Sl 0.009 0.01 
Skin 1000 4.0 0.007 0.01 
Breast 112 79.8 0.139 0.12 
Ovary 11 9.9 0.017 — 
Bladder 43 16.7 0.029 0.04 
Thyroid 33 12.7 0.022 0.04 
Bone marrow 42 61.5 0.107 0.12 
Other solid 144 113.5 0.198 0.12 
Gonads (Heritable) 20 25.4 0.044 0.08 
Total 1715 574 1.000 1.00 


à Included in wr for Gonads. 


b Brain and Salivary glands also each assigned wr = 0.01. 


Wrixon (2008), 


Table 2. Recommended radiation weighting factors. 


Radiation type Radiation weighting factor, wg 


Harrison (2021), 


Table 1. Summary of ICRP publication 103 nominal cancer risks and detriment for uni- 


(from table A.4.1, publication 103, Annex A). 


Photons 


Electrons and muons So, for calculating the absorbed dose, We need 


Protons and charged 


ala is to get the deposited energy in a certain volume 


Neutrons 
2.5 + 
Wa — 504 
2.5 + 


and its mass! 


D = = 


m 


We need to know how determine the 


“CRP 1 
"US NRC 


3 
Š 
0 
u 
o 
5 
£ 
o 
Gİ 
= 


0 
.000001 = .0001 01 1 
Neutron Kinetic Energy (MeV) 


w is still controversial, see, e.g., 


deposited energy 


DUTO 107 TL VI.) 
Other solid 144 113:5 
Gonads (Heritable) 20 25.4 
Total 1715 574 


a Included in wr for Gonads. 


b Brain and Salivary glands also each assigned wr = 0.01. 


9-107 
0.198 
0.044 
1.000 


form whole-body exposure to gamma rays for the whole population, 0—84 years of age 


Tissue weighting 
factor, wr 


0.04 
0.12 
0.12 
0.04 
0.12 
0.01 
0.01 
0.12 


a 


0.04 
0.04 
0.12 
0.12 
0.08 
1.00” 


سے 


1 
2 
3 
7 
5 
6 
1 
8 


Weighting factors w, and w, 


Wrixon (2008), 
doi: 10.1088/0952-4746/28/2/R02 


Photons 
Electrons and 


"d 
رر دا‎ dd 


0 
000001 0001 01 


Neutron Kinetic Energy (MeV 


w. is still controversial, see, e.g.,this presentation 


Harrison (2021), 
0.1088/1361-6498/abe548 


and detriment for uni- 
jon, 0-84 years of age 


Tissue weighting 
factor, wr 


Gonads (Heritable) 20 25.4 
Total 1715 574 


4 Included in wr for Gonads. 
^ Brain and Salivary glands also each assigned wr = 0.01. 


Wrixon (2008), Harrison (2021), Ga! 


Table 2. Recomp 


amm Let's Think on how the total deposited enersy 
Photons 


Ei can be obtained. What do we need to do?‏ ست 


Protons and cl 


ks and detriment for uni- 
ion, 0-84 years of age 


e Tissue weighting 


Alpha particle factor, Wr 


Noia To do that, we need to include methods that T 
allow us to perform "user" actions when: 0.12 


0.12 
0.04 
0.12 
0.01 

^ e GAUserRunAction: the run starts/end 012 
e G4UserEventAction: the event starts/end a 


e G4UserStepAction: the step starts/end 0.04 


17 0.12 


ہے 


/ — 5 
۱۳ =D 5 


9.198 0.12 
Gonads (Heritable) 20 25.4 0.044 0.08 
Total 1715 574 1.000 1.00” 


7 b. : 
600001 6001 0 ٦ a Included in wr for Gonads.  ح»‎ ۰ هم‎ — 
Neutron Kinetic Eneray (MeV) ; b Brain and Salivary glands also each assigned wr = 0.01. 
تا‎ AS Still controversial, SEE, ©ga) 
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G4UserRunAction (and Event and Step) CA. 


Allow tp get control before/after a run. Provides a run-object (G4Run) 
o It is needed when you want to take actions before/after the run 
starts/ends (before first event processing/after last event 
processed) 


Imagine you want to store all the hits in a certain volume (e.g. detector) 


In our case, for the sake of completitude, let’s use a root file instead of an 
ASCII (text) file 


We shall use the g4root standalone libraries for not depending on root installation 


#include "g4root.hh" xo 


Done. Now, all your root knowledge can be included in your app 
G4UserEventAction and G4UserRunAction are similar but acting for events and steps 


o OW FEN دم‎ 


There 


#ifndef RUN HH 
#define RUN HH 


we go.. run.hh/.cc 


#include "run.hh" 


MyRunAction: :MyRunAction() {} 


MyRunAction: :-MyRunAction() (|) 


#include "G4UserRunAction.hh" 


#include "g4root.hh" 


libraries 


class MyRunAction 


public: 


y; 


MyRunAction() 


~MyRunAction( 


virtual void 


//always! root geant4 standalone 


public G4UserRunAction { 


r 


); 


BeginOfRunAction(const G4Run*); 


virtual void 1 


#endif 


EndOfRunAction(const G4Run*) ; 


void MyRunAction::BeginOfRunAction(const G4Run*) { 


G4AnalysisManager *root = G4AnalysisManager :Instance(); 
root-»OpenFile("doses.root"); 

// create the NTuple 
root->CreateNtuple("doses", "doses"); 

// information to be stored in columns 

// root->CreateNtupleIColumn ("fEvent") ; 
root-»CreateNtupleDColumn("fEDep"); 

// root->CreateNtupleDColumn ("fmass"); 

// root->CreateNtupleDColumn ("fAbsDose"); 


root->FinishNtuple(0) ۶ // close the NTuple 


void MyRunAction: :EndOfRunAction(const G4Run*) { 


G4AnalysisManager *root = G4AnalysisManager :Instance() ; 


root->Write(); // always write before to close 


root->CloseFile(); 


O ON © کیپ‎ + EN HL 
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ifndef EVENT HH 
define EVENT HH 


include "G4UserEventAction.hh" 
include "G4Event.hh" 


include "g4root.hh" 


include "run.hh" 


class MyEventAction : public G4UserEventAction { 
jobs ditions 
MyEventAction (MyRunAction*); 
^MyEventAction(); 


virtual void BeginOfEventAction(const G4Event*); 


virtual void EndOfEventAction(const G4Event*); 


void AddEDep (G4double EDep); 


private: 


G4double fEDep; 


#include "event.hh" 
MyEventAction::MyEventAction (MyRunAction*) { 
fEDep = 0.; 


MyEventAction: :~MyEventAction () 1) 


void MyEventAction: :BeginOfEventAction(const G4Event*) { 
// we want to get the deposited energy for each event, then 
// we need to put it to 0. each time the event starts. 
// Otherwise, we will get the total accumulated Ed (later) 
// comment if you want the total accumulated energy 
// first test in an event basis 


fEDep = 0.; 


void MyEventAction: :EndOfEventAction(const G4Event*) { 
G4cout << "Energy deposition: " << fEDep << G4endl; 
G4AnalysisManager *root = G4AnalysisManager::Instance(); 
root-»FillNtupleDColumn(0, fEDep) ; 
root-»AddNtupleRow (0); 


void MyEventAction::AddEDep(G4double EDep) { 
fEDep += EDep; 
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G4UserStepAction (new files: stepping.hh/.cc) (mx 


ifndef STEPPING HH 
define STEPPING HH 


include "G4UserSteppingAction.hh" 
include "G4Step.hh" 


include "construction.hh" 


include "event.hh" 


class MySteppingAction : public G4UserSteppingAction { 
publics 
MySteppingAction (MyEventAction *eventAction) ; 
^MySteppingAction(); 


virtual void UserSteppingAction(const G4Step*); 


private: 

MyEventAction *fEventAction; 
Pe 
#endif 


#include "stepping.hh" 


MySteppingAction::MySteppingAction (MyEventAction *eventAction) { 


fEventAction = eventAction; 


MySteppingAction::-MySteppingAction() 1) 


void MySteppingAction::UserSteppingAction(const G4Step *step) { 
// deposited energy is stored at each step 
G4double stepEDep = step->GetTotalEnergyDeposit(); // this is 
for all volumes 


fEventAction->AddEDep (stepEDep) ; 


Entries 1000 
Mean 9.583 
StdDev 0.9081 


10 


10 


o OM FEN Aa 


We need to know if the particle is in the sensitive vol; 


At ‘construction’ we need to define our sensitive volume, and we need to build a 
function to export this logical volume 
// in the class definition 
private: 
G4LogicalVolume *fSensitiveVolume; 


Py "bg. clus comimos (eina ie meta 


G4LogicalVolume *GetSensitiveVolume() const { return fSensitiveVolume; } 
// define the logical volume in the construction.cc 


fSensitiveVolume = logicMuscle; 
Now, all the information about the current position of the particle is in stepping 
// *volume is the LogicalVolume where the particle is located 


G4LogicalVolume *volume = 


A il me (E‏ ابص تد یں سم دم و محمد 


// we need to get an object including the volumes we constructed 


const MyDetectorConstruction *detectorConstruction = static cast<const 


MyDetectorConstruction*» (G4RunManager: :GetRunManager () ->GetUserDetectorConstruction )( 7 


// and *fSensitiveVolumeNis our selected sensitive volume 


G4LogicalVolume *fSensitiveVolume = detectorConstruction->GetSensitiveVolume )( ۶ 


E, in the whole world 


d 


Entries 1000 
Mean 9.583 
StdDev 0.9081 


E 


d 


in the muscle 


3 4 5 6 


fEDep 
Entries 000 
Mean 


p 
Std Dey/0-9154 


Slightly lower 


10 


fEDep 
Entries 1000 
Mean 0.02217 | 


0 


2 1 -+ 70 
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Success! Now, D=fEDep/sensitiveVolumeMass 


There are several ways to do this. As this is a very beginner course, let's use a 
practical but no so efficient method 
For other ways, check the 'Blexample' included in the Geant4 source directory 


First we need to introduce a new method and attribute at MyEventAction class 


void GetMass (G4double mass); 


G4double fMass; 


Then we need to report the mass to the user eventAction (in stepping.cc) 
fEventAction->GetMass (fMass) 7 


Finally, get the mass at every step (in stepping.cc) (very inefficient!) 
G4double fMass = fSensitiveVolume->GetMass() / kg; // This is the correct way to deal with units 


And just calculate and accumulate the dose at the end of every event (event.c 


G4double fAbsDose = 0.; تہ‎ 
it (ass > 0) 


fAbsDose = (fEDep / joule) / fMass; // This is the correct way to deal with units 


O ON O 0 + EN نم‎ 
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We are done my friends 


1000 e 10 MeV 

Energy deposition: 9560.96 MeV. 
Mass: 0.13125 kg. 

Dose: 1.167 116-00 Gyo — 1127 ۷ 


1000 photons 10 ۷ Ç J | 


Energy deposition: 535.392 MeV 
Mass: 0.13125 kg 
Dose: 6.90615e-10 Gy — 0.690 nGy 


too IA 


RF RF RF PR‏ با 
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Wait! Before you go, remember the ٩ - 


e How can we evaluate the shielding effect? 


e Let's try 


N O <a O UCI + EN P 
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Wait! Before you go, remember the shielding? 


e How can we evaluate the shielding effect? 


e Let's give it a try 


G4Material *lead = 


mişsie—>E imeOrBuileMaterlal (MEM Bo“); 


G4Box *solidShield = new G4Box ("solidShield", 0.30*m, 0.30*m, 0.015*m); 


G4LogicalVolume *logicShield = new G4LogicalVolume (solidShield, lead, "logicShield"); 


G4VPhysicalVolume *physShield = new G4PVPlacement (0, G4ThreeVector(0.*m, 0.*m, 0.0375*m), 


logicShield, "physShield", logicWorld, false, 0, true); 
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remember the shielding, electrons 


1000 e 10 Mev 

no shielding 

Energy deposition: 9560.96 MeV 
Mass: 0.13125 kg. 

Dose: 1.16711e-08 Gy. > 11.7 ۷ 


1000 e 10 Mev 

lead, 3cm 

Energy deposition: 28.70 MeV 
Mass: 0.13125 ke. 

Dose: 3.05302e-11 Gy. > 30.5 ۷ 
Shield effect: -1/400 
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remember the shielding, 


1000 photons 10 ۷ 

no shielding 

Energy deposition: 535.392 MeV 
Mass: 0.13125 kg 

Dose: 6.90615e-10 Gy > 690 pGy 


photons 


1000 photons 10 ۷ 

Lead, 3cm 

Energy deposition: 168.04 MeV 
Mass: 0.13125 kg 

Dose: 2.0513e=10 — 205 pGy 
Shield effect: 1/3.4 


Wait! Before you 


Original dose (no shielding): 690 pGy 


1000 photons 10 ۷ 

Lead, 3 cm 

Energy deposition: 168.042 MeV 
Mass: 0.13125 kg 

Dose: 2.0513e-10 Gy > 205 pGy 


go, remember the shielding? 


1000 photons 10 MeV 1000 photons 10 MeV 

Lead, 5cm Lead, 1cm 

Energy deposition: 62.77 MeV Energy deposition: 382.918 MeV 
Mass: 0.13125 kg Mass: 0.13125 kg 

Dose: 7.66311e-11 Gy > 76.6 pGy Dose: 4.6743e-10 Gy > 467 pGy 
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Conclusive remarks, 1 


e Together we have developed: 
o a simple, not so efficient, but complete Geant4 application 
o it is a powerful tool for performing Geant4 simulations 
o it is a nice template for building your own applications 
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Conclusive remarks, 2 


e This course is based on several Geant4 courses available online. ` 
e My recommendation is to read the all the docs and check at least: 


e https: //www.youtube.com/playlist?list-PLLybsCU6QCGWgzNYOVOSKen9vqg4KXeVL — YouTube 


e https: ithub.com/mnovak42/Geant4-Beginner-Course/tree/master < GitHub 


ns/V10.7/html/index.html — Official guide for G4 apps developers 


I acknowledge Mustafa Schmidt (gphysics matters) and Mihaly Novak for 


their wonderful resources and courses 
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Conclusive remarks, 3 


I am hoping to continue developing my Geant4 course and surely 1 Re l 


introduce changes in the future. 


I will tag this version of the course at GitHub for your future reference 
(Think in a repository as a movie. A tag is picture taken at a certain 
moment of the development) 
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Homework, sorry, you also have to work ;-) 


e Important homework 


o Play a lot with these codes. Try different geometries, 
materials, particle beams, energies, and physics lists. 
Compare the changes and take note of your observations 
Analyze the effect of different shielding materials, thickness 
and positions. 

o Try to build composite shieldings (thin layers stacked of 
different materials) 

o Reproduce the findings on the theoretical courses here: stopping 
power, particle ranges, mass absorption, .. 


e Official homework 
o I will send the final assignation in mattermost later this week 
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Take home messases GA! 


Geant4 is a wonderful toolsuite for simulate the interaction of 
radiation with matter. 


IMHO, G4 analysis tools are not so good, so it I recommend to 
produce root/csv outputs, and analyze them in root/python 
correspondingly 


IMHO, developing a new expertise require: 


1. A simple but functional example . 
۱ this course 
2. Understand the jargon 


3. Google (or Bard/ChatGPT) freely available 


4. Eager to learn it depends on you 
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It was nice to share this course with ۰ (SA 


// “Geant4 for Beginners: a 


crash course” ended here, so 


delete course; 


Thanks for participate! 
Hope to see you soon! 


gasoreyh 


